AXIAL ANOMALY IN THE PRESENCE OF 
THE AHARONOV-BOHM GAUGE FIELD 



Paola Giacconi *, Stephane Ouvry **, Roberto Soldati * 

Q>^ ! * Dipartimento di Fisica "A. Righi" Bologna, sezione I. N.F.N. Bologna, Italia 
^ ** Division de Physique Theorique 

^ ; IPN.F-91406 Orsay, France and 

00 ' LPTPE, Universite Pierre et Marie Curie F- 75252 Paris 

> 

cn 

o 

\o 
o 



ABSTRACT: We investigate on the plane the axial anomaly for euclidean 
Dirac fermions in the presence of a background Aharonov-Bohm gauge potential. 
The non perturbative analysis depends on the self-adjoint extensions of the Dirac 



p ' operator and the result is shown to be influenced by the actual way of under- 
^ . standing the local axial current. The role of the quantum mechanical parameters 
involved in the expression for the axial anomaly is discussed. A derivation of the 
^ ! effective action by means of the stereographic projection is also considered. 
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1. Introduction 

The behaviour of matter fields in the presence of Aharonov-Bohm (AB), or 
anyon, gauge fields has recently attracted a lot of interest owing to its possible 
application to condensed matter physics. Clearly, the origin of the much studied 
bidimensional anyon model [1] is intimately connected with topological AB inter- 
actions. From the field theoretical point of view, the coupling of 2-dimensional 
fermions with such a kind of gauge fields exhibits some interesting features which 
are worthwhile to be carefully investigated. 

In this note we shall deal with the problem of computing the axial anomaly 
[2] , as well as the effective action, for an euclidean 2-dimensional Dirac field inter- 
acting with a background AB potential, whose flux intensity is given by a, with 
— 1 < a < 0. Besides its own theoretical interest, in Ref.[3] the axial anomaly has 
been shown to be related to the second virial coefficient of an anyon gas, thereby 
giving a connection with in principle measurable thermodynamical properties of 
anyonic matter. Moreover it has also been observed that the 1+1 axial anomaly is 
connected with measurable effect in solid state physics [4] . We shall more precisely 
concentrate on the dependence of the axial anomaly on boundary conditions at 
the location of the flux tube, i.e. on possible self-adjoint extensions of the Dirac 
hamiltonian. The result of our analysis will be fully exhibited, on the plane, for 
the special value a = — i, owing to merely technical limitations. As a matter of 
fact, it should be hopefully clear that there is no reason, in principle, to doubt 
that our main statements still hold true for the whole range of a. 

We flrst treat the model in the 2-dimensional plane and later on we shall 
study its compactiflcation by means of the stereographic projection. 

The starting point is the classical euclidean action. 




(1.1) 



where 



ilj) = i-iy,{df, - ieAfj) 



(1.2) 
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the AB gauge potential being 



A 5^ 1 / 2 2 *^ 1/ o\ 

An = —€ai/C)i/ m \/ Xi -\- X2 = —€ui/—n~, 2 ' (-'-•") 

Here, 7i = cti, 72 = <^2, 73 = ^7i72 = <^3, where the sigma's are the Pauh 
matrices. As it is weU known the field strength is 

F^.^le^.F, F=-^5^^\x) , (1.4) 

which corresponds to an infinitely thin solenoid at the origin; if; and ■0 are inde- 
pendent euclidean spinors. 

The basic tool which allows for the non perturbative definition of the axial 
current and anomaly is the complex power of the Dirac operator. In general, com- 
plex powers of pseudo-elliptic invertible operator on compact manifolds whithout 
boundary do indeed exist under very general hypotheses [5]. Since, however, we 
are on the whole 2-dimensional plane and in the presence of field strengths with 
S-\ike singularities, the only way available to set up the complex powers is by 
means of the spectral theorem. To this aim let us consider the eigenvalues and 
eigenfunctions of the Dirac operator . The crucial point is that, in the presence of 
an AB gauge potential, a symmetric Dirac hamiltonian is selected if one considers 
a domain spanned by spinors which are regular at the origin, but then complete- 
ness of the eigenstate basis is lost. In order to find a complete ortho normal basis 
which diagonalizes the Dirac operator it is necessary to consider its self-adjoint 
extensions [6] [7] . 

To this concern let us choose polar coordinates (r, (j)) on the plane; the 
eigenvalue problem becomes 



did did 
(7i[cos(/)— - -sm(f){— + ia)] -F72[sin0— + -cos^(— + ia)])'i/;x{r, 0) 

(1.5) 

where /x is a suitable mass parameter to fix the scale of the eigenvalues. If we 
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rescale the spinor wave functions as 



1 



^(i^)(^),i(n+l) 



(1.6) 



with n e Z, A e R — {0} and ^ = /ir, we get the eigenspinors regular at the origin 
when n ^ 0, namely 



V'A,±n(C,0) =\l^ 



(±z)J±.(|A|e)e±^-'^ 
s(/n(A)J±(,+i)(|A|Oe^(i±")'^ 



(1.7) 



here n e N, being the Bessel function of order j/(±n) = ±n + a. 

On the other hand, the partial waves corresponding to 1^(0) = a can not 
be both regular at the origin unless completeness of the eigenfunctions is lost 
[6]. Then one has to consider the self-adjoint extensions of the Dirac operator 
by means of the standard Von Neumann method of the deficiency indices [8]. 
This leads to a one-parameter family a; G R, whose domain is given by 
V{I^J = + + e''^iP^-'))\i; e C(°°)[0,+oo] nifl([0,oo]),/?e C,u;e R}, 

where 




(1.8) 



Ki, being the Basset-McDonald function, J\f a normalization constant and H2 the 
Sobolev space 

p27v poo 

Hl{[Q,^]) = {^\ / # / ^KiAmiA) < 00} . (1.9) 

Jo Jo 

The corresponding eigenfunctions for u — a can be written in the form 



47r(l + sin^(|A|) cos z/vr) 



X 



icos^Ja(|A|e) -isin^J_«(|A|^) 



sgn{X) cos^J(i+,)(|A|0+sin^J_(i+,)(|A|0 



i\M) 



(1.10) 



where 



tan^(|A|) = |A|2°+itan 



(1.11) 
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We would like to notice that the cigcnfimctions in eq.s (1.7), (1.10) are improper 
eigenfunctions, since they belong to eigenvalues of the continuous spectrum. They 
are suitably normalized according to theory of the distributions, viz. 



hm / ^dU #<(|AlK,0)V'n.(|A2|^,0) = 5n,n,(5(Ai-A2) • (1.12) 

Jo Jo 

Moreover, in order to obtain the correct normalization as in eq. (1.12), one 
has to put the contribution at the origin equal to zero, thereby finding the relation- 
ship of eq. (1.11). It should be stressed that, rather than the requirement on the 
domain X>(^^), which actually involves normalizable states, it is just the condition 
(1.12) of being a complete orthonormal family of improper eigenfunctions, leading 
eventually to eq.s (1.10), (1.11) (for the angular momentum component n = 0). As 
a matter of fact, the property of having a complete orthonormal basis turns out 
to be necessary and sufficient for an operator to be self-adjoint. Furthermore, we 
see that, for any value of cu, the spectrum is purely continuous and is provided by 
the whole real line, but the point A = 0, the zero modes being absent since they 
are not orthonormalizable. 

2. The £Lxial anomaly on the plane 

Once the eigenvalue problem has been solved, we are able to set up the 
complex power by means of the spectral theorem. The complex power of the 

< Cl,0l|d6,02 >= i^-s(/a;;6,01,6>02) 



dimensionless operator I^^ = [ ) is defined by the kernel 



+ 



OO r /.CXD , 

2 J] / dA A->A,n(ei, <l>im n(6, 02) + {u ^ -n) 
n=l 

/•oo , 

2 / dA A-«^i';)(6,0i)^i"o^n6,</'2) , 

^0 



(2.1) 



which can be analytically extended to a meromorphic function of the complex 
variable s; the key property is that the kernel of the complex power is regular at 
s = 0. In particular, the value of its trace over spinor indices, on the diagonal 



iCii 4^1) = (^2, 02)7 can be explicitely evaluated either in the case uj = 0,71, — 1 < 
a < or in the case a = w G R. Actually, to the aim of computing the 

anomaly, it is more useful to write down the explicit form of the traces of the axial 
kernels, namely 

tr[^3Ks {Iu,=oy, e, e)] = ^-s (a, e) + (-«, C) 

^^.s-2 r(^)r(a + i-f) ^ (2.2) 

+ — f L^i—L — 1 £i \-(a\ — > a + 1) 

where we have set 

s iao- ^^cs-^ r(^)r(a + 2-f) 

^-^"'^^-2.^ 0F(.-i)(.-2)r(f)r(a + f) ' ^'-'^ 

whereas 

tr[73i^.(/.;e,0]la=-i = 



r(f) V r(^) 

Now we are ready to obtain the local forms of the axial anomaly, as it 
arises from the non perturbative definition of the fermionic axial current associated 
to the invertible operators I^^. As a matter of fact we notice that, since the 
kernel K_s{I^; ^,^) of the complex power is a well defined tempered distribution, 
depending meromorphically upon the complex variable s, one can properly define 
the euclidean averages of the vector and axial currents, respectively, by means of 
point-splitting as well as analytic continuation [9], namely 

< >= e < tr[7^V;(a;)V;t(x)] >= 



lim lim e tr[^aK-s{Iuj'i x,x + e)] 

s— >1 e— >0 

+00 /.+00 



(2.5a) 



lim hm J2 / e *K7mV'S(^)V'S + ^)] 



where < • > means euclidean average and 

< >= e < tr[7p73^(a;)V^'''(x)] >= 

lim lim e tr['ya'y3K_s(I^; x,x+ e)] = 

+ 00 p + OG 



(2.56) 



lini lim / dX e tr['y^-f3ii;l^l{x)ip!^l\x + e)] 



From the above definitions of tlie averaged local currents it is straightforward to 
show, taking eq. (1.5) into account, the quantum balance equations, namely 



<d^j^^\x)>^Q , (2.6a) 

testing the gauge invariance of the definition in eq. (2.5a), whereas 

< >= 2ie lim tr[73 Ks{I^; x, x)] = A^'^\x) (2.66) 

leads to the definition of the local axial anomaly, once the topology has been 
chosen in taking the limit s — > 0; we shall discuss below this delicate matter. 
First we consider the <S'-topology, namely we study 



lim (Px tr[y3 Ks{Iuj;x,x)] f{x) , 
s-^o J 



for any rapidly decreasing function / e 5(R^). 

If we take the limit s — > in the sense of the distributions, from eq. (2.2) 

we get 

' (fx A^°'''\x)fix) = -a , (2.7) 



where / is a suitable test function belonging to »S(R^) normalized to /(O) = 
we notice that the above result, in full agreement with the one of Ref.[3], actually 
corresponds to the usual result, viz. A^'^''^\x) = — ^e^j^F^j^(x) as a distribution. 

On the other hand, we can start from eq. (2.4), with a = and try to 
take the limit s ^ 0, for any u), in the <S'-topology. Let us consider indeed in the 
RHS of eq. (2.4) the term proportional to sincu. Then from the identity 



5' -lim (e---^-ILsi-\^)) = 

s^O \ S J 



(2.8) 

1 d djlnp^ 

2^de^ dC ^ ' 

it immediately follows that the 5'-limit does not exist unless a; = 0, tt and, in this 
case, eq. (2.7) still holds with a = —\. The above analysis strongly suggests that 
the definition of the anomaly, as the limit in the sense of distributions of the axial 
kernel, leads to select u; = 0,7r as the correct physical choices for the self-adjoint 



extensions of the Dirac operator. The same conclusion has been claimed in the 
literature [7] from a quite different point of view, namely by demanding that one 
does not have any additional contact interaction at the origin beyond the point- 
like magnetic field. Nevertheless it is worthwhile to mention that some dual point 
of view leads to a different outcome, when we consider the limit s — > in the 
natural topology of R — {0} and only afterwards continue the result to »S'(R^). 
As a matter of fact, a non vanishing result is obtained in this case for a; 7^ 0,7r 
and, when a = — |, we can compute explicitely 

iG sm tu 

liuiliuietr[^3K_s{Iu,-,x,x + €)]\^^_i = r^O . (2.9) 
As a consequence, there is a unique continuation in »S'(R^) which reads {a = — |) 

noo /»27r roo /'27r 

cdcj^ dci>m,c/>)A^^\o- 1 cdcj^ ^'^ISpj ' (2.10) 

where / is a test function belonging to »S(R^). We recall the definition of the 
tempered distribution 

^ ^ :{d! + di){lnrf + C{fi)5^^\x) , (2.11) 



[e] ~ 2/.2 

C{iJ,) being an arbitrary function of the rcgularization mass parameter. It should 
be stressed that the freedom (actually up to In//) within C(/x) amounts to the 
natural requirement 

' ' ' ; (2.12) 



consequently we can think about C{p) in terms of a scaling function, in the sense 
that, when fx pfx, we have C(/i) C{pi^) = C'(a*) + 27rlnp. 

Now some remarks are in place to comment this result. On the one hand, 
we observe that eq. (2.10) involves two further quantum mechanical arbitrary 
parameters, u) and the scaling factor C{iJ,), beyond the flux intensity a = — | of 
the classical background infinitely thin solenoid. In particular, if the test function 
vanishes at the origin, where the field strength is concentrated, still a nonvanishing 
contribution survives, of a purely quantum mechanical nature, which depends 
upon the parameter of the self-adjoint extensions. On the other hand, we recall 



that in the case of smooth gauge fields the local axial anomaly turns out to be 
proportional to the field strength itself, up to the presence of zero modes (see Refs. 
[3] [9] [10] and the last part of the present paper). In the presence of the singular 
AB gauge potential, the local axial anomaly still exhibits the above feature if we 
define the axial current as a limit in the »S'-topology, and, moreover, the self- 
adjoint extension is fixed. Alternatively, the local axial anomaly is different from 
zero even outside the infinitely thin solenoid, and the parameter of the self-adjoint 
extensions is free, if we define the axial current and anomaly as in eq.s (2. 9), (2. 10). 
We could say that this latter quite interesting feature is closely reminescent of the 
AB eff'ect. 

In conclusion some physical, a priori measurable quantity [3] [4], such as 
the local axial anomaly, in the present singular case, appears to be reproduced by 
the standard Schwinger's [11] form A{x) = — ^e^i/F^i/(a;), the parameter of the 
self-adjoint extensions being fixed to cu = 0, tt, if the regularized axial current is 
defined as a limit in the topology of tempered distributions. Nonetheless, from the 
mathematical point of view, some inequivalent construction of the singular axial 
current might be considered, which eventually leads to a non standard form of 
the axial anomaly and to the freedom in the choice of the self-adjoint extensions. 
Strictly speaking, the above statements holds true in the special cases we have 
explicitely worked out. Nevertheless it is highly presumable that the same features 
still appear in general, although the explicit proof does not seem to be presently 
available. The physical content of this last approach will be discussed elsewhere. 

3. The effective action 

Another important issue is the calculation of the effective action. To this 
purpose it should be noticed that the standard gauge invariant Schwinger formula 
on the plane [11] does not make sense, owing to the singular nature of the AB 
gauge potential. Nevertheless, a meaningful definition of the effective action may 
be obtained [9] [12], through the transition to the compact case by means of the 
stereographic map from the plane to the open 2-sphere (without the north-pole) , 
the origin of the plane coinciding with the south-pole. The open 2-sphere is 
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embedded in R"^ in such a way that a convenient choice of the coordinates reads 



y cos (p 
ysiiicf) 

l + y^ 

2 
3 



Xi = 2a 
X2 = 2a 
X3 = a 

1 - 

X^X^ + Xl = a'', y 



(3.1) 



r 

a 



It is immediate to obtain the form of the zwei-beine, namely 



1 + y 



2 "^CM ' 



Kiy) = ^-^Sa, 



(3.2) 



//'27r poo 



(3.3) 



The invariant measure becomes 

Ty2)2 

Within the above conformal coordinate system Xf^, the eigenvalue problem for the 
covariant Dirac operator takes the form 



' f 1 



^"^^ V'A(a^) + i-V'A(a^) = , 
a 



(3.4) 



where 



V^ = 5^-ieA^(ar) . (3.5) 

If we pass to polar coordinates and perform a conformal mapping on the spinors, 
namely if we set 



I > 



(3.6) 



the eigenvalues problem becomes equivalent to the system of coupled differential 
equations 



y 



y 



dy ""^^ "z(l + y2)^A 
ay i[l + y^) 



(3.7a) 
(3.7&) 
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where i^(±n) = ±n + a, z/(0) = a. In order to solve eqs. (3.7) it turns out to be 
useful to perform the change of variables 



j, I dKp, cf>) = j^J d<t> £^ (3.8) 



^ — 2 

y 

Hence, the eigenvalues problem takes the form 

il-p')^-l^x[''^+^Xl^/l^x[''^ = , (3.9a) 

, (R) 

{l-p')^ + {i^ + l)x'i''^+i\iV^^x'i''^=0 , (3.96) 



with Xi ^ and 



(l-p^)^-.xi^)=0 , (3.10a) 



{i-p^)^ + {v + i)xf^ ' ^^-^^^^ 

for the zero modes. 

We would like to notice that the covariant Dirac operator on the open 2- 
sphere, up to the conformal mapping (3.6), turns out to be essentially self-adjoint 
with domain given by square integrable functions on the interval p e [—1, 1], the 
derivatives at p = ±1 being understood in the sense of the distributions. In 
order to fulfil the above requirement, only the value n = is allowed, a feature 
which leads to the removal of the degeneracy with respect to the orbital angular 
momentum. The complete set of orthonormal eigenfunctions reads 

a 

P. 



^' ^^''^^^ 27^1,1^) ^1 ' ^3.11a) 



l + g 



X<«'(p.*) = f^(l-P)(i±^) ' P/_T'"°'(P) , (3.116) 

with eigenvalues Xi = I = 1,2, .... The sign indetermination is related to 
the global euclidean axial symmetry if; i — > e^^'^i/j. Here Pi'^'^\p) stands for the 
normalized Jacobi polynomials [13], namely 




p{a,i3)f . / /!(2/ + a + (3 + l)r(/ + a + (3+l) (g,/?) . . 

KP)-\\ 2"+/^+ir(/ + a+l)r(/ + /3 + l) ' ' ^3^2) 

tll!(i - p)-"(i + p)-^^[(i - pr^\i + pf^' . 
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The noteworthy feature is the appearence of zero modes, namely 

xi^\p) = -^Vsm{-Tva) ' (3.13a) 

xi''\p) = , (3.136) 

which indicates the non trivial topological behaviour of the AB potential. By 
the way, the reason for the vanishing of the right-zero-mode comes from the 
requirement of orthonormality for the complete set of the eigenfunctions of the 
self-adjoint Dirac operator on the open 2-sphere. 

Let us come to the explicit evaluation of the quantum effective action and 
of the corresponding axial anomaly. As it is known, the 2-dimensional classical 
action for the Dirac spinor field on the open 2-sphere is conformally equivalent to 
the corresponding quantity on the complex plane, namely [12] 

Wpiane = Wsphere + / d^X (jj{x)Au{x) , 

. ^2 (3-14) 

where W stands for the logarithm of the regularized determinant of the Dirac 
operator. As a consequence, taking eq. (3.6) into account, we can define the regu- 
larized determinant of the flat Dirac operator by means of the zheta-regularization 
of the associated self-adjoint operator of eqs. (3.4), which has a discrete spectrum. 
Actually, as the spectrum includes the null eigenvalue, we have to set [9] 

oo 

C(my;5)^5](Z2)-^ = C(2.) , (3.15) 
1=1 

where 

. (3.16) 

and consequently 

Wpiane = -'^T^det' {iaf) = -^C{iaf,s)\s=o = ~ln27r ■ (3.17) 

whereas 

W^spftere = ^(l-ln2)- ^ln27r . (3.18) 
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The local axial anomaly is now arising from the variation (see Refs. [8] [14]) 

/ d^x A{x)f{x) = 2ie [ d^x ^-^^ -ex^d^fix) , (3.19) 
J J dAi^{x) e 

and leads, nota bene, to a one-parameter family of local expressions, the parameter 
being given by the compactification radius a, namely 

oo 

Aix) = 2ielim V(/2)-^ trMi{-)i;ji-)] . (3.20) 

s ^■U ^ ^ CI CI 

Unfortunately, it does not seem that the RHS of the last equation could be set 
into a close analytical form, just as in the general case of the continuous spectrum 
on the plane. It should be stressed that the standard Seeley-DeWitt asymptotic 
expansion can not be safely used to handle the RHS of eq. (3.17), since the origi- 
nal theory is strictly defined on the open 2-sphere without poles; in particular the 
fermion eigenfunctions are singular (not meromorphically) at the poles. Neverthe- 
less, one can try to approximate the present singular case by means of a smooth 
"vortex-like" gauge potential, namely one can consider 

Of X 

A^{x)= 5'-lini -e^^^- — ^- — ^ , (3.21) 

^ cr^o e ^ xi + X2 + cr^ 

the limit being taken at the end of the anomaly calculation (actually this corre- 
sponds to a regularization of the singular statistical AB interaction) . Within this 
framework, the anomaly can be computed from the standard methods [9], [10] and, 
in the decompactification limit a oo, it reproduces the usual Schwinger's result, 
the topological zero modes being disappeared. It should be noticed that this result 
is in agreement with what we have previously found on the plane, if the anomaly 
is coherently understood as a limit in the topology of tempered distributions (see 
eq. (2.7)). In this sense we argue that the results on the plane, we have previously 
discussed, should be of general character. 

In conclusion we find that the effective action can be computed exactly (see 
eq. (3.17)) starting from the stereographic map, whereas the local axial anomaly 
can be explicitely evaluated on the plane, where the effective action does not exist 
in the present singular case. Here we have presented the detailed calculation of 

13 



the local axial anomaly on the plane only for particular values of the parameters. 
Obviously, a technical effort should be attempted in order to work out exactly the 
quite general cases for the effective action as well as for the anomaly. It would 
also be very interesting to investigate other kinds of physical situations such as, 
in particular, the occurrence of confining potentials and/or smooth boundaries, to 
further understand their influence on the axial anomaly and the effective action. 
Moreover we think to analyze the eventual relevance of eq. (2.10) in the evaluation 
of many body hamiltonian, virial coefficients and other physical effects. 
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